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We simulate the terahertz emission from laterally-biased InGaAs and InP using a three- 
dimensional carrier dynamics model in order to optimise the semiconductor material. Incident 
pump-pulse parameters of current Ti:Sapphire and Er:fibre lasers are chosen, and the simulation 
models the semiconductor's bandstructure using parabolic F, L and X valleys, and heavy holes. The 
emitted terahertz radiation is propagated within the semiconductor and into free space using a model 
based on the Drude-Lorentz dielectric function. As the InGaAs alloy approaches InAs an increase 
in the emitted power is observed, and this is attributed to a greater electron mobility. Additionally, 
low-temperature grown and ion-implanted InGaAs are modelled using a finite carrier trapping time. 
At sub-picosecond trapping times the terahertz bandwidth is found to increase significantly at the 
cost of a reduced emission power. 

PACS numbers: 42.72.Ai, 73.20.Mf, 78.20.Bh, 78.47.+P 



I. INTRODUCTION 

The ability to create and detect single-cycle pulses of co- 
herent terahertz (THz) radiation has enabled diverse stud- 
ies in condensed matter physics. Recent examples include a 
study of how surface plasmon polaritons propagate on semi- 
conductor gratings,. determining the low-energy vibrational 
modes of oligomers,'^ and measuring the time- dependent inter- 
subband absorption in semiconductor quantum wells. Such 
studies have become possible using high power, broadband 
THz radiation emitters based on the ultrafast separation of 
photoexcited carriers in a photoconductive switchi- In this 
communication we investigate how the choice of semiconduc- 
tor alters THz emission, by employing a carrier dynamics sim- 
ulation to model InGaAs and InP photoconductive switches. 

The use of a low-bandgap semiconductor (such as 
Ini_xGaxAs with an alloy fraction x < 0.47) would enable 
cheaper, more stable, turnkey THz spectroscopy setups based 
on erbium fibre (Er:fibre) lasers, which can produce pulses as 
short as 65 fs at a central wavelength A — 1.55 /imi^ Combin- 
ing a compact system of this type with communications wave- 
length optical fibres would enable THz pulses to be generated 
and detected in small, extreme, or otherwise inaccessible en- 
vironments (e.g. in an endoscope or in the cavity of a pulsed 
magnet) . 

While GaAs-based photoconductive switch technology is 
well developed-, there are fewer InP and InGaAs-based THz 
emitters described in the literaturej^*^*^*— *— . Baker et. ali~ 
reported THz emission up to 3 THz from a low-temperature 
(LT) grown Ino.3Gao.7As photoconductive switch, using a Nd 
glass laser. Additionally, Suzuki and Tonouch;^" have mea- 
sured the THz emission from photoconductive switches made 
from Ino.53Gao.47As implanted by 340 keV, 1 x 10^^ cm~^ Fe 
ions. They used an Er:fibre laser with central wavelength 
1.56 /xm and pulse duration 300 fs, and measured THz power 
up to ~ 2 THz. Finally, Mangeney et. ali— performed a simi- 
lar experiment on Ino.53Gao.47As implanted by Au^ ions, with 
a carrier lifetime of 0.9 ps. 



In order to investigate suitable semiconductors for photo- 
conductive switch THz emitters it is beneficial to adopt a 
modelling approach. The effect of altering one parameter (or 
more) of the semiconductor material or laser can be investi- 
gated rapidly, and an optimum combination sought. While 
analytical models of THz emission are relatively straightfor- 
ward to interpret, they often require phenomological assump- 
tions (e.g. about carrier mobilities). In this paper we there- 
fore employ a three-dimensional semi-classical carrier dynam- 
ics simulation to model Ini_xGaxAs and InP photoconductive 
switch THz emitters. 



II. SIMULATION DETAILS 

The carrier dynamics model used herein is an extension of 
one used previously to model the THz emission from InAs 
and GaAs surfaces under a magnetic fieldr" ion-implanted 
GaAs surfaces— and GaAs photoconductive switches^ A 
set of 10^ pseudoparticles comprising extrinsic and photo- 
generated carriers and fixed ions are used to simulate the 
semiconductor. At each 5fs step in time the model numeri- 
cally finds the three-dimensional potential due to the charge 
density, subject to the appropriate boundary conditions for 
the surface and Schottky photoconductive switch contacts. 
Carrier- carrier, carrier-phonon and carrier-impurity scatter- 
ing mechanisms are included^ The simulated particles are 
within a box of size xxyxz = 6x6x4: /im^ , subdivided into 
a grid of 64 x 64 x 32. The contacts of the photoconductive 
switch are at —3 < a; < — 1 nm and 1 < a; < 3 fim. A con- 
stant bias voltage of 5 V between electrodes was used through- 
out this work, which corresponds to typical experimental field 
strengths. The simulation starts at a time t = — 0.7ps, to al- 
low the extrinsic carriers to equilibrate before the arrival of 
the incident pulse, which has peak intensity at time t — 0.0 ps, 
and is centred at a;, y = 0. 

Although the simulation assumes a parabolic band struc- 
ture for simplicity, in the modelling of THz emission from 
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narrow band-gap semiconductor surfaces the inclusion of the 
L-valley was found to create a more realistic distribution of 
carrier energies^ In photoconductive switches the average 
carrier energy is greater than for surface emitters, and inter- 
valley (e.g. r ^ L) scattering by absorption or emission of 
an optical phonon becomes more significant. Further, in the 
range 0.4 < a; < 0.7 the X-valley energy gap of Ini_xGaxAs 
falls below that of the L-valley. The simulation was therefore 
extended to include the X valley in addition to the F and L 
valleys. 



A. Properties of InGaAs 

Values of semiconductor parameters for Ini_xGaxAs at 
a particular alloy fraction x were generated automatically 
within the simulation from the bowing parametes— When 
no bowing parameter was found in the literature, linear inter- 
polation between InAs and GaAs was used. The two-mode 
TO phonon energies for InGaAs of Ref»-- were employed. For 
consistency the intervalley deformation potentials were taken 
from Ref^, the only reference found with values for GaAs, 
InAs and InP. 

An inspection of the properties of InGaAs suggests poten- 
tial advantages and disadvantages as a THz emitter when 
compared with GaAs. Firstly, a lower F-valley effective mass 
should result in a higher electron mobility, and greater THz 
power. Secondly, the intervalley deformation potential Dr-L 
is greater for InAs than GaAs^-, resulting in a larger inter- 
valley scattering rate. When electrons scatter to the L or X 
valley they lose kinetic energy, contributing to the decelera- 
tion of carriers that produces the negative peak in Ethz- The 
lower absorption coefficient of InGaAs at 1.55 fim compared 
with GaAs at 800 nm (0.8, 1.2 x 10*^ m~^) is a disadvantage, 
because a lower a results in a lower bandwidth for defect-laden 
semiconductors, owing to the photogeneration of fewer carri- 
ers within the damaged depth. In addition, with a smaller 
a the emitted THz radiation will on average travel further 
through the semiconductor (when the refiection geometry'' is 
used), and will be more significantly affected by absorption 
and dispersion. 

Finally, since Ini_xGaxAs exhibits two-mode TO phonon 
behaviour^^ the THz emission spectra should exhibit strong 
absorption peaks near the TO phonon energies (e.g. for 
x=0.47 the two modes have frequencies of 6.8 and 7.6 THz). 
This will limit the spectroscopic use of Ini_xGaxAs photocon- 
ductive emitters in the frequency range near and between the 
TO phonon modes of InAs and GaAs (6.4 and 8.1 THz). The 
following section describes how the two-mode TO phonon- 
polariton dispersion of InGaAs was included in the simula- 
tion. 



TO phonon modes is given by: 

eTo(a;,x) = e^ -I- —-^ " ■ , 2 



l-c.Va.2^,-ic.r/a.2^i ^ ^ 

where ej, are the low and high frequency dielectric con- 
stants at alloy fraction x. The angular frequency of the GaAs- 
like (InAs-like) TO phonon is tJx^i (i^x=o)- The damping 
rate was assumed to be F = 0.06 ps~^ for both InAs-like and 
GaAs-like TO phonon modesf2S In order to estimate the thick- 
ness of semiconductor through which the emitted THz must 
propagate we calculated the maximum THz electric field from 
each 'slice' parallel to the surface. Taking the mean of this 
distribution results in a weighted propagation thickness of 
~ 0.5 /im. Assuming that the emitted THz radiation is col- 
lected in the refiection geometry (to minimise absorption and 
dispersion in the semiconductor material^) 0.5 fim is a reason- 
able estimate of the propagation distance of THz radiation 
in an InGaAs photoconductive switch. We therefore propa- 
gate i?THz through this thickness of semiconductor, and into 
free space, resulting in a characteristic reduction in spectral 
amplitude of emitted THz radiation close to a TO phonon 
mode, and an enhancement just above the frequency of an 
LO phonon (as seen in Figs. Illjt . 



C. Laser parameters 

The simulation uses incident pulse parameters characteris- 
tic of current ultra-short mode-locked lasers (i.e. a Gaussian 
spatial and temporal shape, and a transform-limited Gaus- 
sian energy distribution) . For the case of a 10 fs Ti:sapphire 
laser (with central wavelength A = 800 nm, AA = 80 nm, 
typical beam power Pcxp ~ 400 mW and a repetition rate 
R = 75 MHz) we use a simulation power of 10 fiW and a 
Gaussian spot of standard deviation a^^y = 0.5 ^m in or- 
der to obtain the same photon fiux as experimentally achiev- 
able. While amplified Er:fibre lasers have been demonstrated 
with pulse durations of 65 fs (A = 1550 nm, AA > 100 nm, 
Pcxp = llOmW, R = 67MHz),^ sub-30fs pulses widely tun- 
able over the wavelength range 1130 to 1950 nm may be gen- 
erated by coupling such Er:fibre laser light into a highly non- 
linear fibrei— For excitation around 1.55 /im we use P = 6 jj.W 
and (Ta;,j, = 0.5 fim, chosen to give the same photon fiux as for 
the Ti:sapphire, so that a comparison between Ino.53Gao.47As 
switches excited by 1.55 /^m radiation, and GaAs excited by 
800 nm can be made in Section ITTTl With the aim of providing 
a quick reference to the laser parameters used we refer below 
to the 10 fs, A = 800 nm Ti:sapphire just described as laser A, 
and the 65 fs A = 1.55 ^im Er:fibre as laser B. 



B. THz propagation within photoconductive switch 

The radiated THz electric field Etuz was calculated in the 
far-field approximation from the simulated current density J 
using the relation Ethz oc (1 -|- ^/ej'^dJ /dt^ where the di- 
electric function e was modelled using the well-known Drude- 
Lorentz formula for TO phonon-polaritonsf-- which for two 



III. SIMULATION RESULTS 

In Fig.Qwe compare the simulated Ethz for InP and GaAs 
photoconductive switches excited by 10 fs pulses from laser 
A with that obtained for InAs excited by 65 fs pulses from 
laser B. These parameters were chosen to enable a comparison 
between THz systems based on current state-of-the-art lasers. 
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In order to enable a comparison between the results presented 
we normalize Ethz in the time and frequency domains to 
the peak simulated values for a GaAs photoconductive switch 
excited by laser A (with carrier trapping time Tc — 100 ps) , 
denoted Eq and Eq respectively. Similarly the emitted THz 
power (defined as the square of ETUz{t) integrated over all 
times) is normalized to Po, the power for the GaAs emitter. 

The peak value of Ethz for InAs is larger than for GaAs, 
despite the Er:fibre laser's power being lower than that of 
the Ti:sapphire (the photon flux is the same for both lasers). 
This can be attributed in part to a higher fraction of incident 
photons having above-bandgap energy. As the inset of Fig. Q 
indicates, almost 100% of photons from laser B have energy 
greater than i5r(InAs) = 0.35 eV, while the above-bandgap 
fraction for laser A is smaller (~80%). The fact that InAs 
has a lower effective mass in the F- valley than GaAs also con- 
tributes to the power increase. Qualitatively, a smaller effec- 
tive mass produces a larger mobility fi = er/m* , which should 
result in a larger transient current and therefore power. How- 
ever, this argument neglects changes in the carrier scattering 
time r, and we return to this point in Section fill Gl 

Turning now to -Ethz for InP, we observe a slightly greater 
peak time-domain and spectral amplitude than for GaAs, de- 
spite the larger (lower mobility) of InP. This is again 
related to the energy distribution of the incident infrared 
pulses. Defining the excess energy of the incident pulse as 
A = E^ - Er we see that A = 0.21 eV for InP, while 
A = 0.13 eV for GaAs. A greater percentage of incident pho- 
tons are therefore absorbed in InP (~91%) than GaAs (80%), 
contributing to a larger Ethz- 

The simulated Ethz for InAs exhibits an extra oscillatory 
component after t = 0.30 ps, due to the onset of charge- 
density oscillations that can radiate THz radiation. These 
plasma-type oscillations can be observed directly in the simu- 
lation as a change in electron density between time steps (Fig. 
^c)). At t = 0.00 ps the rapid increase in electron density is 
due to photoexcitation, and produces the first positive peak 
in Ethz- By t — 0.10 ps there is an increased electron den- 
sity at the anode, while the first negative peak in Etuz is at 
t = 0.2 ps (Fig. ^a)), and is primarily a result of the rapid 
momentum scattering of electrons. The subsequent change in 
Ethz for t > 0.3 ps can be identified as resulting from circu- 
lar regions of altering charge density (Fig. [H^c). t — 0.55 ps). 
Plasma oscillations have been observed previously in the sim- 
ulation of THz emission from InAs surfaces, and become 
less noticeable as x — > 1 due to GaAs's lower mobility. 



A. Carrier trapping time 

Higher bandwidth terahertz emission can be achieved by in- 
troducing defects into a semiconductor using a technique such 
as high energy ion implantation, or low-temperature (LT) 
semiconductor growth. In both cases, carriers are trapped 
on sub-picosecond time scales, with carrier trapping times Tc 
as short as 550 fs in LT-Ino.3Gao.7As (Refi-) and 300 fs in 
Ino.53Gao.47As implanted with 2MeV Fe ionsi— When Tc be- 
comes comparable with the duration of a single-cycle THz 
pulse (tc < 1 ps) bandwidth improvements can be observed. 
This effect may be modelled using an exponential reduction 
in the photoexcited carrier density, with a time constant Tc. 



X (|im) 
-1 1 




Frequency v (THz) 



FIG. 1: (Colour online) (a) Simulated Etuz for GaAs (dashed 
line) and InP (dash-dotted line, offset in time by -1-0.2 ps 
for clarity) photoconductive switches excited by laser A, and 
InAs (solid line) excited by laser B. Both InP and InAs 
have a larger peak Sthz than GaAs. (b) Fourier transform 
of Etuz (i) from (a), for InAs (solid) and GaAs (dashed 
line). The spectrum for InP (not shown) is substantially 
the same as that for GaAs (with its TO phonon absorp- 
tion at 9.2 THz). (Inset) Schematic energy distributions of 
a transform-limited 65 fs pulse at A =1.55 /^m (0.8 eV) and a 
10 fs pulse at A =800 nm (1.55 eV). Vertical lines at 0.73 eV, 
1.34 eV, 1.42 eV and 1.71 eV respectively indicate Er for 
Ino.53Gao.47As, InP and GaAs, and El for GaAs. (c) Change 
in electron density between simulation time steps in the plane 
of the surface {x ~ y), averaged over the depth z (see text for 
discussion). The black lines at a; = — 1/im and x = 1 fim 
mark the positions of the cathode and anode. 



In Fig. 121 we present the simulated THz electric field Etuz 
and calculated spectra for Ino.53Gao.47As excited by 65 fs 
pulses of 1.55 fim centre wavelength light, as a function of car- 
rier trapping time. At shorter values of Tc the negative peak in 
Ethz arrives earlier, creating an increase in the spectral fuU- 
width-at-half- maximum (FWHM). This trend can be seen in 
Fig.|2Ic): when Tc < Ips there is a rapid increase in FWHM 
at the cost of a reduced power, where the power is defined as 
£'THz(i) summed over all simulation times. In the following 
we chose a value of Tc = 0.3 ps, achievable experimentallyi— 



B. Varying alloy fraction at A = 1.55 fim 

With the intention of optimizing the semiconductor for 
excitation at 1.55 fim we modelled the THz emission from 
Ini_xGaxAs at alloy fractions x in the range < x < 0.7. 
The results of this set of simulations are shown in Fig.|3 The 
emitted power increases towards InAs to over four times larger 
than that of GaAs excited by laser A, while the FWHM re- 
mains constant. When x > 0.45 both the power and FWHM 
drop dramatically, since the excess incident pulse energy A 
becomes small and then negative, and fewer photoexcited car- 
riers are created. 
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FIG. 2: (Colour online) (a) Simulated -Ethz for Ino.53Gao.47As 
with Tc = 100 ps (solid line) and Tc = 0.3 ps (dashed line), us- 
ing 65 fs pulses with laser B's parameters, (b) Fourier trans- 
form of (a). In (c) the power (circles) and FWHM (crosses) 
are plotted for Tc = 0.1 — 100 ps, showing a significant re- 
duction in power and a bandwidth increase as Tc decreases 
below 1 ps. The power is normalized to that of GaAs with 
Tc — 100 ps (excited by laser A), Po- The dotted vertical line 
indicates Tc =0.3 ps, used in subsequent figures. 
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FIG. 3: (Colour online) (a) and (b) Simulated Ethz for InAs 
(solid lines) and Ino.53Gao.47As (dashed lines) with Tc =0.3 ps 
and excited by 65 fs, A = 1.55/im pulses from laser B. (c) 
Power (circles) and FWHM (crosses) of simulated Ethz as a 
function of x. The dotted vertical line indicates x=0.47. 



C. Varying alloy fraction at constant excess pulse 
energy A 

Simulations were run over the full range of alloy fraction 
X with A = 0.07 eV and constant photoexcited carrier den- 
sity. Pulsed Er:fibre lasers widely tunable in the range 1130 
to 1950 nm (1.1 to 0.64 eV) have been recently developed;— 
which should allow experimental studies of this kind over a 
significant range of alloy fraction (0.37 < x < 0.77). 

As Fig. Ilja) shows, the emitted THz power increases to- 
wards x=0, as found in Section llll Bl By extracting the to- 
tal carrier momentum scattering time r from the simulation 
(summed over all scattering mechanisms, averaged over all 
electrons in the F valley, and taken a.t t — 0.00 ps) we estimate 
the high-field electron mobility ^(x) ~ er(x)/mp(x). The rel- 
ative increase in ^i(x) towards InAs is comparable with the 
enhancement in experimental low- field mobility (Fig. Elb)), 
and results from a decrease in both the carrier-carrier and 
carrier-LO phonon scattering rates. 

Shortly after the peak intensity of the incident pulse the 
electric field applied across a photoconductive switch is com- 
pletely screened by the build-up of the THz-emitting dipole. 
This occurs in a characteristic screening time ts ~ 100 fs for a 
semi-insulating GaAs photoconductive switchi By calculat- 
ing the electric field at the incident spot position {x,y — /im) 
as a function of time, it is possible to extract the screening 
time ts , defined as the time at which the electric field becomes 
zero. Figure 01b) indicates that ts increases as the Ga alloy 
fraction x 1, as a result of the lower electron mobility. It 
should be noted that a shorter screening time does not neces- 
sarily mean a reduction in the THz power emitted: a higher 
mobility can produce both a lower ts and a greater power. 

This simulation set was repeated for incident pulse dura- 
tions of 10 fs and 100 fs. In the 10 fs case the total power 
emitted (Fig. 0Ja)) is smaller than for 65 fs across the entire 
range of x, as fewer carriers are photoexcited. However, for 
100 fs the power is again lower than the 65 fs case, because 
when the later carriers are photoexcited the electric field is 
already partly screened. How the THz power and bandwidth 
alter as a function of incident pulse duration was discussed 
for semi-insulating GaAs in Ref.— . The FWHM of the 100 fs 
spectra are lower than those of the 10 fs and 65 fs simulations, 
owing to the slow rise in Etuz- Somewhat counterintuitively, 
the 10 fs spectra do not have greater FWHM than the 65 fs 
over the whole range of x. This can be explained by noting 
that the FWHM measures the width of the spectrum, and not 
the high-frequency spectral tail. Thus while greater spectral 
amplitude is obtained at high frequencies for 10 fs excitation 
(e.g. Fig. 0^b)), the amplitude at low- frequency (which dom- 
inates the FWHM) is larger at incident pulse durations of 
65 fs. 



While these results should be relatively straightforward to 
verify experimentally, the density and kinetic energy of pho- 
toexcited carriers alter significantly when x is varied and A 
kept constant. To try to assess better how the material choice 
alters THz emission (independent of the laser's parameters) 
we performed a further set of simulations at constant excess 
excitation energy A, as reported in the next section. 



IV. CONCLUSION 

In summary, we have simulated the terahertz emission from 
Ini_xGaxAs photoconductive switches over the entire alloy 
fraction range < x < 1. Two methods for increasing the 
THz power from photoconductive switches were identified. 
Firstly, a larger electron mobility resulted in a larger THz 
power as the Ga alloy fraction was decreased. This sug- 
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FIG. 4: (Colour online) (a) Power (solid lines) and 
FWHM (dashed lines) of simulated THz electric field from 
Ini_xGaxAs (t^ — 0.3 ps), for incident pulse durations of 10 fs 
(crosses), 65 fs (triangles), and 100 fs (squares). The photon 
energy above the bandgap and the incident photon density 
were kept constant as the alloy fraction x was varied, (b) The 
average F valley electron mobility (solid line) extracted from 
the simulation (see text for description), can be compared 
with the (low carrier energy) Hall mobility (dashed line). The 
time ts in which the applied electric field becomes screened 
(squares) increases towards GaAs. 



gests that semiconductors with high mobility /low effective 
mass may be used to improve THz power, for example InAs 
(mp — 0.022me) or InSb (mf = 0.014me). Experimentally, 
a high resistivity or semi-insulating InAs may be required to 
produce a photoconductive switch with a low dark current. 
Secondly, increased THz power can be obtained by judicious 
choice of incident pulse energy, so that the majority of pho- 
tons are absorbed. 

The spectra calculated from the simulation indicate that 
InAs photoconductive switches have greater power than GaAs 
over the frequency range — 15 THz. At higher frequen- 
cies, the initial rise in Etuz (dominated by the incident pulse 
duration) governs the amplitude. This result seems to sug- 
gest that an In As/65 fs Er:fibre laser setup should produce 
THz pulses that are at least as powerful as those from 10 fs 
Ti:sapphire/GaAs setups. 
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